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Using a quantum mechanical three-step model we present numerical calculations on the high- 
harmonic generation from four polyatomic molecules. Ethylene (C2H4) serves as an example where 
orbital symmetry directly affects the harmonic yield. We treat the case of methane (CH4) to 
address the high-harmonic generation resulting from a molecule with degenerate orbitals. To this 
end we illustrate how the single orbital contributions show up in the total high-harmonic signal. This 
example illustrates the importance of adding coherently amplitude contributions from the individual 
degenerate orbitals. Finally, we study the high-harmonic generation from propane (C3H8) and 
butane (C4H10). These two molecules, being extended and far from spherical in structure, produce 
harmonics with non-trivial orientational dependencies. In particular, propane can be oriented so 
that very high-frequency harmonics are favorized, and thus the molecule contains prospects for the 
generation of UV attosecond pulses. 
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I. INTRODUCTION 

In the process of high-harmonic generation (HHG) 
high-frequency coherent radiation is emitted from an 
atom or molecule when exposed to intense femtosecond 
laser light. This process has been the subject of exten- 
sive studies during the last couple of decades due to pos- 
sible applications, including, e.g., generation of coherent 
UV attosecond pulses p], [Jy] , carrier-envelope phase re- 
trieval H and tomographic reconstruction of molecular 
orbitals @, H, 0| • 

Whereas early HHG studies focused on atomic sys- 
tems, it was not until the new millennium that exper- 
imental work on molecules gathered pace. Studies on 
HHG from molecules is motivated by the expectation 
that the more degrees of freedom and non-spherical sym- 
metry as compared to atoms may lead to richer physics 
and a higher degree of control. Hitherto the molecules 
considered have mostly been diatomic or linear systems 
and it is only just recently that experimental results 
on HHG from more complicated molecules have become 
available [8[ . Here we apply a quantum mechanical three- 
step model to investigate theoretically HHG from poly- 
atomic molecules of current interest to experimental- 
ists jgl, [H . In our theory, we use a detailed description of 
the molecular orbitals, obtained from Hartree-Fock cal- 
culations. We illustrate the effect of orbital symmetry 
on HHG by studying the orientation dependence of the 
harmonic signal from ethylene (C2H4). The methane 
molecule (CH4) is used to demonstrate how degenerate 
HOMOs influence the harmonic spectrum. In particular, 
we show the importance of including coherently the con- 
tribution from every single HOMO when calculating the 
harmonic yield from an oriented molecule. We demon- 
strate how the extended and non-spherical molecules 
propane (C3H8) and butane (C4H10) give rise to har- 
monic spectra with a rather complex orientational depen- 



dence. In the case of propane harmonics near the cutoff 
can be selected by orienting the molecule, and since we 
expect such harmonics to be synchronically emitted Q 
orientation of propane is identified as a tool for generat- 
ing attosecond pulses. It is the first time, to our knowl- 
edge, an investigation of the detailed orientational depen- 
dence of the single molecule high-harmonic signal from 
such complex systems is presented. 

The paper is organized as follows. In Sec. [Ill we derive 
formulas for the single-particle respons from a system 
with degenerate HOMOs, and we review the quantum 
mechanical three-step model used to calculate the har- 
monic yield [10|. In Sec. [1111 we briefly discuss some nu- 
merical details for our calculations. Section HVl contains 
numerical results describing the harmonic yield from each 
of the aforementioned molecules. Especially, we treat 
several effects related to the orientation of such systems. 
Finally, we give a summary and the conclusions in Sec.lVl 
[We use atomic units (e = h = m e = ciq = 1) through- 
out.! 



A. 



II. THEORY 



Harmonic yield from a statistical mixture of 
molecules 



The complex amplitude for the emission of harmonics 
with frequency cjhhg polarized along the linear polariza- 
tion vector e is obtained from the Fourier transform of 
the dipole acceleration 



A e (^hhg) = e 



/ 



(1) 



where (d)(t) is the expectation value of the dipole oper- 
ator d of the molecule. The corresponding power density 



2 



is given by [11|, [12[ : 

^e(^HHG) OC |A e (cJ H HG)| 2 . 



(2) 



The quantum state will be a mixed state, since several 
unobserved variables appear in a measurement of the har- 
monic signal. First, typically several electrons from a 
given molecule will contribute significantly to the HHG 
because of degeneracy of the highest occupied molecu- 
lar orbital (HOMO). An experiment recording the to- 
tal harmonic yield from an ensemble of ground state 
molecules cannot distinguish contributions from the in- 
dividual degenerate HOMOs. Second, since perfect ori- 
entation of molecules is unrealistic, the measured har- 
monic yield arises from an ensemble of molecules with 
different orientations. In this section we treat these un- 
resolved degrees of freedom, using and elaborating the 
results of recent work [13[. Initially, the molecule is in 
a stationary thermal state at temperature T, and the 
system is completely characterized by the density ma- 
trix po = exp (—H/kBT)/Z, with partition function Z = 
Tr[exp(—H/kBT)], H the field- free molecular Hamilto- 
nian and ks Boltzmann's constant. We want to resolve 
the molecular initial state on energy eigenstates. The 
molecules we consider are, however, so complex that sev- 
eral simplifying assumptions have to be made in order to 
obtain a practical theoretical formulation. According to 
the Born-Oppenheimer approximation, we can separate 
the electron and nuclear motion, and we will assume only 
the electronic ground state be populated. Next we adapt 
the single-active-electron (SAE) approximation, and in- 
troduce an index A in order to be able to discriminate be- 
tween the degenerate HOMOs of this active electron. The 
nuclear motion consists of rotation and vibration. We 
shall assume that only the vibrational groundstate is oc- 
cupied, whereas the rotation needs to be treated in more 
detail to take account of an oriented molecule. In do- 
ing so, we neglect the rovibrational interaction such that 
rotation and vibration can be treated separately. The ro- 
tation is characterize d by the asymmetric top quantum 
numbers J, r and M [14j . We may, accordingly, specify 
the energy eigenstates by their electronic and rotational 
degrees of freedom, |A) <g> \JrM). When the system in- 
teracts with a laser, the energy eigenstates will evolve 
according to a unitary operator that describes any num- 
ber of orienting pump pulses followed by a short, intense 
laser pulse that drives the HHG [spontaneous decay pro- 
cesses can be neglected on the timescales we consider]. 
The orientation pulses are not strong enough to affect 
the electronic motion appreciably jl5[, while the nuclear 
dynamic can be considered as frozen during the short 
high-harmonic generating pulse [13[. Thus, if the delay 
between the final orienting pulse and the driving pulse is 
denoted by td, we can split the time evolution operator 
according to U(t) ~ ?7hhg(£) ® [/orient (td) f° r a descrip- 
tion of the evolution when the HHG is produced, where 
[JhhgW propagates the electronic part of the molecular 
state during the pulse of the driving laser and [/orient (^d) 
accounts for the evolution of the rotational state of the 



molecule. The energy eigenstates then evolve as follows 



■ JtM (*)> - (U WG (t)\\)) ® (U OTient (t d )\JrM)) 



\Mt)$JTM(td))- 



(3) 



We return to the time evolution of the system in Sec. Ill Bl 
Our current goal is to calculate the expectation value 
of the dipole operator as this enters Eqs. (pQ) and ([2]). 
This evaluation is most conveniently done by expanding 
the energy eigenstates in the position basis |r, 0, 0, x) m 
which the dipole operator is diagonal 



(d(t)} = Tr \p(t)d] = Tr \u(t)p^U ] (t)d 
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dr\^ x (r,(t),e,x^)\ 2 d, 



(5) 



J2 Pjt\$Jtm(^0, X ,U)\\ (6) 

J,r,M 



and Pj T = ex.p(—Ej T /kBT) the Boltzmann weight of the 
asymmetric top wavefunction (0, 0, x\ [^orient (td)\JrM) = 
$jymO, 0, x; td) of energy E Jr . 

Consequently, using Eqs. (pQ) and ([2]), the harmonic 
signal is given by 



Se(^HHG) cx 



E 



2n p7r rlir 
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X G td (0, 0, X)A* (^HHG , 0, 0, X) 



(7) 



with 



(8) 

In general, the degenerate HOMOs will interfere due to 
the coherent sum over A in Eq. ([7|). Nevertheless, there 
are special cases where the degeneracy enters simply as 
a factor multiplying the signal from a single HOMO. For 
example we mention that if the degeneracy is due to spin 
multiplicity, A^s, the sum over different HOMOs yields a 
factor A/| in the signal. Another instance occurs when 
the HOMOs differ simply by a rotation. Then, in the 
case of a randomly oriented ensemble of molecules, i.e., 
Gt d — l/(87r 2 ), it is obvious that HOMOs must each 
give rise to the same complex number when the single 
HOMO amplitude, A^(ujylylGi 0? x)> is averaged over 
all orientations. It then follows from Eq. © that the 
degeneracy will again enter as a factor multiplying the 
signal from a single HOMO. 
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B. Model 

As described in the previous section the system starts 
out in a statistically mixed state composed of the energy 
eigenstates. Each |A) <g> \JrM) state evolves according to 
Eq. (j3]) . In this section we focus on the time evolution of 
the system. 

First, propagating the asymmetric top energy eigen- 
states, | JVM), to obtain the \<&j r M(td))'s at the time td 
of the high-harmonic probing pulse is a numerically de- 
manding task. The issue has already been addressed in 
several studies [HI, [H, [I?}, and in the present work, we 
will assume either Gt d to be simply uniform (no preferred 
orientation) or use an idealized orientational distribution 
to be specified in Sec. HVl 

Next, we focus on the electronic part of the time evolu- 
tion. In the field-free initial state the HOMO wave func- 
tion is conveniently expressed in a spherical expansion in 
the body fixed (BF) frame 

^{r)=Y J FUr)Y l m {r). (9) 

Asymptotically this expression must follow the Coulomb 
form 

V>f» -Y, C t m r Z/K - 1 eM-^r)Yr(r) (10) 

l,m 

with k, = v // 27p, Ip the ionization potential, Z the 
net charge of the molecule when the HOMO electron 
is removed and where the Cj A m 's are fitting coefficients. 
More detail on how to obtain wavefunctions in Eqs. ([9]) 
and ([T0|) is given in Sec. IIII1 We wish to carry out calcu- 
lations in a laboratory fixed (LF) system defined by the 



laser polarization, and for a molecule of arbitrary ori- 
entation LF and BF coordinate axes do not in general 
coincide. Hence, we rotate the BF wave function to ob- 
tain the LF wave function by application of the rotation 
operator 

V^ F (r, 4>, 0, X ) = Dtf, 9, X )^ F (r), (11) 

where the rotation is given by the Euler angles </>, and 
X- Following the conventions of [l4| is the angle be- 
tween the BF z-axis and the LF Z-axis, <p denotes a ro- 
tation around the Z-axis, and finally x denotes a rotation 
around the z-axis. Note that the effect of D is readily 
evaluated in the spherical harmonic basis used in Eqs. ([9]) 
and CLOD 0. 

We consider the case where the driving laser pulse con- 
tains several cycles such that a Floquet approach is suit- 
able. Hence using the Coulomb gauge and the dipole 
approximation, a laser field with frequency uo and pe- 
riod T = 27i/uo is described by the vector potential 
A(t) = Aocos(cut). According to the quantum me- 
chanical three-step model described in [18[ for atoms 
and [To[ for molecules the electronic time evolution, given 
by Uu¥LG(t) [see Eq. (|3j)], consists of a HOMO electron 
being transfered to the continuum via above threshold 
ionization (ATI), i.e., by absorbing a number of photons 
from the driving laser. The electron then propagates in 
the laser-dressed continuum and is eventually, due to the 
periodicity of the laser field, driven back to a recombina- 
tion with the molecule, where it returns to the HOMO. 
Within this model the complex amplitude for the emis- 
sion of harmonics polarized along the unit vector e with 
frequency cjjjhg = Nuj (N integer) is 



h,h m f 2 ,m' 1 m 2 ,mi k C(k) 

Here D l m , m .(</>, x) w ^h i = 1, 2 is the Wigner rotation function [l4|, while 



C(k)> 
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q=K k +A(t) 



(14) 



along with their Wigner rotation functions, are inter- k photons during the ATI-step. In Eqs. ([T3j) and ([T4j) q 
preted as ATI and propagation-recombination ampli- 
tudes, respectively, of a HOMO electron having absorbed 
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and q' are electron momenta and 



S(t) = kuot 



Kk • — - sm(uut) 

00 



^ sin(2ut) (15) 
2oo 



is the quasi-classical action. The index C(k) in Eqs. ([12]) - 
(fT4|) denotes the saddle-points. For each /c the saddle- 
points t'c{k) are defined by the condition S'(t' c ^) = 0, 
and we use the ones with < Re(t f c ^) < T along with 
Im(t' c(/c) ) > 0. The factor (±l) h in Eq. (13| corresponds 
to the limits ±m of the size of the electron momen- 
tum at the saddle-points. The factor 1/Lq(£, ^c(fc)) = 
aa^sincc;^^ — sinc^t) in Eq. (Jlj), with a = ±1 to as- 
sure Re(Lo) > 0, describes the decrease of the ampli- 
tude of the electron wave as it propagates in the field- 
dressed continuum. Also, is the part of the contin- 
uum electron momentum arising from absorption of k 
laser photons during ATI, thus = y/2{ku — I p — U p ) 
with Up = the ponderomotive potential and e^ fc = 
ae A • Finally, in Eq. ([T4j) the function Ffe m2 (q) is the 
radial part of the momentum space HOMO wave func- 
tion, obtained by taking the Fourier transform of Eq. (jTTJ) 
(see [10] for further details). 



to understand the harmonic yield from a whole gas of 
molecules. 

In all results presented below, we calculate the signal 
of harmonics polarized along a linearly polarized 800 nm, 
1.8 x 10 14 W/cm 2 driving laser. As the light is linearly 
polarized the results are independent of 0, the rotation 
around the polarization vector. We consider molecules 
that are either randomly oriented or have been one- or 
three-dimensionally oriented. The orientational distri- 
butions G td [see Eqs. ©-(0)] corresponding to random, 
one-dimensional or three-dimensional orientation are as 
follows 



^random 



1 

1 5(9' -9) 
47r 2 sin 9' 
1 6(9' -9) 



2?r sin 9' 



s(x' - X)- 



(16) 
(17) 
(18) 



Here the BF z-axis is oriented at angle with respect to 
the LF Z-axis in the cases of one- and three-dimensional 
orientation and the molecule is rotated an angle \ around 
the BF z-axis in the case of three-dimensional orienta- 
tion. 



III. CALCULATIONAL DETAILS 

We have determined the Hartee-Fock wave functions 
for ethylene (C2H4), methane (CH4), propane (C3H8) 
and butane (C4H10) in a spherical basis along with 
the asymptotic coefficients using the technique described 
in [l9j]. In Table U we list the ionization potentials I p and 
the asymptotic C z A m -coefficients entering Eq. ([TO]) . Note 
that the results in the case of ethylene differ from those 
in [19[ because our choice of the BF axes in the present 
work follows the convention of [16[ with the xz-plane co- 
inciding with the molecular plane, whereas in pjl the 
molecular plane was chosen to coincide with the ?/z-plane. 



IV. RESULTS AND DISCUSSION 

In this section we present results on the high-harmonic 
yield from ethylene, methane, propane and butane. Since 
we do not include effects of propagation [22, [23[ our cal- 
culations cannot be directly compared with the experi- 
mental results of Ref. [8]. If our main purpose was the 
optimization of the harmonic yield with the object of gen- 
erating attosecond pulses then, surely, phase-matching 
should be taken into account. Nevertheless, in this the- 
oretical work aiming at isolating and illustrating clearly 
the effects of symmetry, degenerate orbitals and orienta- 
tion, we find it reasonable to disregard phase-matching: 
first, because propagation effects can be reduced exper- 
imentally by using a gas jet which is short compared to 
the coherence length [24j |. and second, an understand- 
ing of the single molecular response is needed in order 



A. Ethylene: Effects of the orbital symmetry 

We first present results on the HHG from ehtylene 
(C2H4). This molecule has a non-degenerate HOMO 
which makes the influence of the HOMO on the har- 
monic signal relatively transparent. Additionally, ethy- 
lene is interesting from the point of view that field-free 
one-dimensional alignment has been carried out experi- 
mentally [16] and field-free three-dimensional alignment 
has been explored theoretically [17[. Consequently, the 
theoretical results presented below may, in principle, be 
subject to experimental investigations. 

Figure [1] shows HHG spectra from ethylene at different 
orientations corresponding to the orientational distribu- 
tions of Eqs. ([T6]) and ([T7]) . The overall effect of orienting 
the molecule is a scaling of the spectrum. The reason 
for this scaling is that all harmonics have similar orien- 
tational dependence. We also note the absence of even 
harmonics in the spectrum, which is easily explained from 
the inversion symmetry of the HOMO of ethylene: The 
inversion symmetry means that the HOMO is composed 
of angular momentum states separated by even multiples 
of h (see Table H]). At the same time absorption and emis- 
sion of laser photons change the size of the electronic or- 
bital angular momentum by ±ft, so in order both to start 
off and end up in one of the angular momentum states of 
the HOMO an even number of dipole transitions is nec- 
essary. Hence, the HHG process requires the absorption 
of an odd number of photons followed by the emission of 
a single odd harmonic. 

We have investigated HHG from C2H4 that has been 
fixed in both 0- and %-angles corresponding to the three- 
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TABLE I: The molecular properties of the alkalenes used in this work for evaluation of HHG. I p is the experimental adiabatic 
ionization potential [20]. The remaining numbers in the table give the values of the asymptotic coefficients, C* m , entering 
Eq. Jl0} based on GAMESS calculation [2l| using a triple zeta valence basis set with diffuse sp shells. The three degenerate 
HOMOs in CH 4 are denoted by HOMOl, HOM02 and HOM03. 
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FIG. 1: (Color online) The orientational dependence of the 
harmonic spectrum from ethylene. The absence of even har- 
monics is explained by the inversion symmetry of the HOMO 
(see text). 



of ethylene in the figure. The meaning of the Euler an- 
gles was explained below Eq. ([9]). In the model used to 
simulate HHG an electron has to escape along the laser 
polarization axis [cf. Eq. ([T3])]. This is impossible, if 
the polarization axis lies along the nodal plane, which is 
the reason for the vanishing harmonic signal, when ei- 
ther = 0° (180°) or x = 0° (180°, 360°). When the 
molecule is rotated the nodal plane is removed from the 
polarization axis of the laser and the strength of the har- 
monics increases. As seen from Figs. [2^b) and (c) the 
harmonics peak at different values of the Euler angle 0. 
The varying positions of the peaks arise from different 
orientational behavior of the ionization and propagation- 
recombination steps making up the HHG process [cf. the 
discussion below Eq. ([T2 |) -(p^ )) ]: As the electron escapes 
along the polarization direction the ionization is maximal 
when lies in between 0° and 90° [l9[. The propagation- 
recombination step, however, is optimized when = 90°, 
but the width of the peak depends on the harmonic order. 
These observations account for the different orientational 
behavior of the harmonics shown on Figs.[2jb) and (c). 

We note, in passing, that a set of data as the ones 
presented in Figs. E^b) and (c) for a full range of har- 
monic energies would, in principle, allow a tomographic 
reconstruction of the HOMO [5] . 



dimensional orientation given by Eq. ([T8]) . Figure [2] 
shows representative results of the calculations. In or- 
der to understand the results, we also show the HOMO 



B. Methane: Interference of degenerate HOMOs 

We now turn to the harmonic yield from the methane 
molecule (CH4). Methane has three degenerate HOMOs 
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FIG. 2: (Color online) The orientational dependencies of the 
harmonics reflect the HOMO of ethylene. In particular the 
harmonics vanish, when the laser polarization coincides with 
the nodal plane of the HOMO. Panel (a) shows the geom- 
etry of ethylene (C2H4) along with an isocontour for the 
HOMO. Dark shading (red online) denotes a negative sign 
of the HOMO wave function. Light shading (brownish on- 
line) denotes a positive sign. Directions of the BF axes are 
shown. We always choose the center of mass as the origin 
of the BF coordinate system. Panels (b) and (c) present the 
dependencies of the 21st (H21) and 29th (H29) harmonics on 
orientation as given by Euler angles and %. 
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FIG. 3: (Color online) The geometry of the methane molecule 
(CH4) along with isocontours of the three degenerate HO- 
MOs. The signs of the HOMO wave functions are indicated 
by the coloring, where dark shading (red online) denotes a 
negative sign and a light shading (brownish online) denotes a 
positive sign. The coordinate system shows the directions of 
the BF axes. Note that the HOMOs differ from one another 
only by a simple rotation. 



as shown on Fig. [3l and we can use this molecule to 
demonstrate the effect discussed below Eq. (|7J), i.e., the 
absence of interference effects from randomly oriented 
molecules, when degenerate HOMOs differ only by a ro- 
tation. To this end we have compared the total harmonic 
yield with the yield from a single HOMO when the orien- 
tational distribution is as given by Eq. (fTB|) and confirmed 
that the results agree except from an overall scaling fac- 
tor. This is illustrated on Fig Ufa). 

Next, we have carried out calculations of HHG from 
methane that has been one-dimensionally oriented with 
the BF z-axis at some fixed angle relative to the po- 
larization direction. The orientational distribution used 
is given by Eq. (|T7|). We do not show the harmonic 
spectrum in this case, since it does not differ much in 
structure from Fig. Ufa). This is probably due to the 
fact that methane is rather small and compact, which 
makes it spherical-like after x-averaging. Consequently, 
no structure is revealed by the electrons, not even the 
most energetic, and the harmonics exhibit the same over- 
all ^-dependence. 

On Fig. 0|(b) we show this typical angular dependence 
of the harmonics. In the figure the upper curve shows the 
signal when the coherence between the individual HO- 
MOs is correctly accounted for [cf. Eq. ([7j)]. The other 
curves in the figure show the unphysical signals from 
each HOMO. Clearly, this figure illustrates that there is a 
strong interference between the single HOMO amplitudes 
in the angle resolved signal. We may understand the sin- 
gle HOMO signals in Fig. ffl(b) from the structure of the 
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FIG. 4: (Color online) The figure illustrates the interference 
of high- harmonics coming from different orbit als as prescribed 
by Eq. J7|. The physical system used is methane (CH4). 
Panel (a) shows the total harmonic yield from randomly ori- 
ented methane. It is just a multiple of the harmonic yield 
from a single HOMO. Even harmonics vanish as a result of 
the orientational averaging. Panel (b) illustrates the detailed 
^-dependence of the 11th harmonic (Hll). Although not mea- 
surable, we also show the single HOMO signals. The calcu- 
lated signals have all been normalized to the maximum total 
yield of the 11th harmonic. The figure underlines the impor- 
tance of adding coherently the amplitude contributions from 
each HOMO, A = 1,2,3. 



HOMOs: First, we explain the dips. The vanishing sig- 
nals of HOMOl and HOM02 at = 0° are explained by 
the fact that in this case the polarization vector points 
along the BF z-axis, hence coincinding with the nodal 
planes of these HOMOs, as seen from Fig. [3l Further- 
more, the electron causing HHG must escape along the 
polarization axis to Eq. (jT3j) , and we conclude that these 
two HOMOs cannot generate harmonics for = 0°. At 
6 = 90° the polarization axis is directed into the BF xy- 
plane which is the nodal plane of HOM03 (see Fig. [3]) 
and consequently HHG from HOM03 is excluded. Sec- 
ond, we remark that the single HOMO-yield at a given 
value of 6 is obtained by averaging amplitudes from all 
degrees of rotation around the BF z [cf.-axis © and (|I7|)]. 
It is therefore obvious from symmetry that the results of 
HOMOl and HOM02 must be identical. 



C. Propane and butane: Effects of orientation 

In the following we consider propane (C3H8) and bu- 
tane (C4H10). We begin with numerical results on 
the propane molecule. Figure [5] shows the HOMO of 
propane along with the harmonic signal in randomly 
oriented and one-dimensionally oriented scenarios cor- 
responding to orientational distributions from Eqs. (jT6j) 
and (|T7|). The dependence of the HHG spectrum on ori- 
entation amounts to more than just a simple scaling, e.g., 
the 19th harmonic is alternately above and below the 
neighboring odd harmonics depending on the orientation 
of the propane molecule. We ascribe this to the fact 
that propane is a rather extended and open-structured 
molecule. In general, when the BF z-axis is close to the 
LF Z-axis the HHG is suppressed, because the almost 
vanishing wave function along the polarization direction 
[see Fig. [5(a)] makes the ATI amplitudes in Eq. ([T2]) 
small. Notice also that the even harmonics disappear 
at an orientation of 90°. This is due to the fact that 
every photon absorption at this orientation changes the 
projection of the electron angular momentum on the Z- 
axis with ±ft, and it is seen from Table [J that an odd 
number of photon absorptions is necessary in order to 
start off and end up in the HOMO. 

We remark that harmonics of order 15 to 27 are sup- 
pressed for = 90°, which means that alignment is a tool 
for favorising the harmonics close to and above the semi- 
classical cutoff energy at I p + 3.17U P — 29 x Huj. These 
harmonics are known to be well phase-synchronized in 
the case of diatomic molecules [9[, and this may also 
hold true in the case of propane. Furthermore, phase- 
synchronized harmonics are synchronically emitted and 
the superposition of such harmonics constitute the basis 
for the generation of attosecond pulses. 

Figure [6] illustrates the HHG from butane. As with 
propane we observe a complicated behavior of the indi- 
vidual harmonics upon the orientation. Opposite to the 
case of propane, there is no particular orientation that 
favors harmonics near the semiclassical cutoff energy. 
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FIG. 5: (Color online) The figure clearly illustrates how the 
open and extended structure of propane leads to a harmonic 
spectrum with a complicated orientational dependence. Panel 
(a) shows the geometry of the propane molecule (C3H8) along 
with an isocontour of the HOMO. Dark shading (red online) 
denotes a negative sign of the HOMO wave function and light 
shading (brownish online) denotes a positive sign. The xyz- 
axes of the BF frame is shown. Panel (b) illustrates the de- 
pendence of the harmonic spectrum on the orientation of the 
propane molecule. At the orientation 6 — 90° even harmonics 
are suppressed as discussed in the text. 



V. CONCLUSIONS AND OUTLOOK 

In the present work we have discussed general issues 
related to the effects of molecular symmetry, degeneracy 
and orientation in HHG. To this end, we have consid- 
ered the high-harmonic signal from several polyatomic 
molecules, namely ethylene, methane, propane and bu- 
tane. In the case of ethylene we have shown how the de- 
pendence of orientation reflects the HOMO: The strength 
of the harmonics is increased when the laser polariza- 
tion is directed away from the nodal plane, but the de- 
tailed orientational dependence differ from one harmonic 
to another. We have used methane as an example of a 
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FIG. 6: (Color online) Panel (a) shows the geometry of the 
butane molecule (C4H10) and an isocontour of the HOMO. 
The sign of the HOMO wave function is indicated by the col- 
oring, where dark shading (red online) denotes a negative sign 
and light shading (brownish online) denotes a positive sign. 
We also show the coordinate system of the BF frame. Panel 
(b) shows the harmonic spectrum corresponding to different 
orientations. As in the case of ethylene the HOMO is inversion 
symmetric which excludes the presence of even harmonics. 



molecule with degenerate HOMOs. Then harmonic am- 
plitudes of different HOMOs need to be added coherently, 
and interference effects are in general unavoidable. Al- 
though, some information about the individual HOMOs 
can be extracted by orienting the molecule we illustrated 
the importance of including all HOMOs when calculat- 
ing the harmonic signal from one-dimensionally oriented 
methane. Finally, the propane and butane molecules 
served as illustrations of HHG from extended structures, 
and the individual harmonics carry their own character- 
istic orientational dependence. In the case of propane 
harmonics near the cutoff can be promoted by means of 
orientation and a selection of these energetic harmonics 
is of interest to attosecond pulse generation. 

While we expect the conclusions drawn above to be 
fairly model independent we would like to mention a few 
open theoretical questions that should be addressed in 
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the near future. In the present paper a relatively sim- 
ple model using several approximations, including the 
strong-field approximation and a saddle-point approxi- 
mation, has been applied to describe HHG. The model 
is formulated in length gauge, but at this level of ap- 
proximation it is still to be settled wether another gauge 
would be more appropriate [25|, [26[. In addition, a fu- 
ture challenge is to take into account the finite duration 
of the harmonic driving laserpulse. This is necessary in 



order to investigate further the possibility of using large 
molecules as sources of attosecond pulses. 
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